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We investigate theoretically the shielding of the electromagnetic field fluctuations by s-wave and 
d-wave superconductors within the framework of macroscopic quantum electrodynamics. The spin 
flip lifetime is evaluated above a niobium and a bismuth strontium calcium copper oxide (BSCCO) 
surface, and the screening effect is studied as a function of the thickness of the superconducting 
layer. Further, we study the different temperature dependence of the atomic spin relaxation above 
the two superconductors. 
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I. INTRODUCTION 

Recent advances in the trapping of neutral atoms near 
microstructured surfaces P, [l] , colloquially referred to as 
atom chips, has proved these structures to be promis- 
ing for the coherent manipulation of cold atom clouds. 
Applications of atom chips include the imaging of the 
electromagnetic field near metallic and dielectric surfaces 
d, 0], the accurate measurement of gravity [1, @j and 
atom interferometry 0, @ . One of the central challenges 
proves to be the interaction between the trapped atoms 
and the hot surfaces in their vicinity, which leads to 
thermally-induced spin relaxation dynamics and attrac- 
tive Casimir-Polder forces. The atomic lifetime is limited 
by the transitions between trapped and untrapped states 
induced by fluctuating electromagnetic fields arising in 
conducting materials ^, [l^ . 

The origin of the electromagnetic field fluctuations lies 
in the finite resistivity (and the finite temperature) of 
conducting materials as a consequence of the fluctuation- 
dissipation theorem. Superconductors have a (almost) 
negligible resistivity and their thermal fluctuations are 
attenuated with respect to a dielectric or a metal. An 
important enhancement of the atomic lifetime is expected 
with the adoption of superconducting surfaces, and both 
theoretical and experimental investigations have recently 
turned their attention to superconducting materials [111 - 

In this article, we investigate the shielding of the near- 
field noise in the vicinity of a dielectric substrate provided 
by two different superconducting films. Casimir-Polder 
potential and thermally induced spin flips are among the 
effects which may be altered when a superconducting 
layer is placed above a dielectric structure, and we will 
focus on the latter. 

Previous theoretical investigations on superconducting 
atom chips have been concerned with the increase of the 
spin flip lifetime and vortex detection via spin dynam- 
ics Three different descriptions of superconduc- 
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tivity have been adopted to estimate the spin flip rate 
[l8l |: the two-fluid model, the Bardeen-Cooper-Schrieffer 
(BCS) theory and the Eliashberg theory. The Eliashberg 
theory provides the most appropriate model for the elec- 
tromagnetic energy dissipation in superconducting ma- 
terials, whereas the other two theories ignore the strong 
modification of either the imaginary part or the real part 
of the optical conductivity in the superconducting state. 
Nevertheless, the spin fiip lifetime evaluated with the 
two-fluid model gives remarkably accurate estimates of 
the Eliashberg results [l^ . 

Quantum mechanical approaches for spin flip lifetime 
calculations above superconducting atom chips have con- 
sidered only isotropic conventional (s-wave) supercon- 
ductors. In this paper we introduce a comparison be- 
tween the screening properties of a conventional s-wave 
superconductor such as niobium, and an unconventional 
d-wave superconductor such as bismuth strontium cal- 
cium copper oxide (BSCCO). The main difference be- 
tween these two classes of materials is the way in which 
the electrons interact with each other to form Cooper 
pairs p^ . In conventional s-wave superconductors, the 
two electrons of a Cooper pair are in a state with zero 
total spin and total angular momentum, such that their 
state is isotropic and the superconducting phase is ho- 
mogeneous with an uniform penetration depth A. 

In d-wave superconductors, an unconventional pairing 
takes place and the state of the paired electrons has a zero 
total spin while the total angular momentum is 2h (hence 
the name). Such superconducting state is anisotropic 
giving rise to layered superconductors, which are charac- 
terised by an anisotropic permittivity. This results from 
the in-plane penetration depth A|| being different from 
the longitudinal penetration depth Aj^. 

As an example, we report the values of the penetra- 
tion depth for the two superconducting materials consid- 
ered in this work. For zero temperature, niobium has 
a penetration depth of A(0) — 35 nm [23|, whereas the 
transverse and in plane penetration depth of BSCCO are 
A^(0) = 100 ^m and A||(0) = 300 nm, respectively \2l\. 

Furthermore, d-wave superconductors have rather high 
transition temperatures. For example, for BSCCO one 
finds transition temperatures of up to Tc = 90 K while 
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for thin niobium films one has Tc = 8.3 K (compared 
to Tc = 9.3 K for bulk niobium). The origin lies in 
the crystal structure of unconventional superconductors 
[19'] . The interlayers act as a reservoir of charge carriers 
that combine into Cooper pairs and, depending on the 
concentrations of oxygen atoms into the material, the 
transition temperature increases from small values up to 
a maximum value of 90 K. 

Not only the value of the transition temperature, 
but also the temperature dependence of the penetration 
depth differs in both types. In conventional supercon- 
ductors, the deviation of the penetration depth A from 
its value at T = can be approximated by the empirical 
formula [191 
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while in d-wave superconductors, A||(T)/A||(0) increases 
with T as [li, Hi 
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In present experiments on superconducting atom chips, 
the superconductor is deposited on an insulator, as its 
contact with a metal can perturb the superconducting 
properties. However, in order to emphasize the screen- 
ing properties of a superconducting medium, we consider 
the superconducting film to be ideally placed above a 
metal substrate. We emphasize here that all the results 
obtained are valid only for the superconductors in the 
Meissner state. The Meissner state is observed when a 
magnetic field in any point of a superconductor is be- 
low the first critical field, which is 140 mT for niobium 
and 13 mT for BSCCO at 4.2 K. In an higher field, the 
magnetic field penetrates into the superconductor (type- 
II) in the form of vortices. In this so-called mixed state 
a superconductor exhibits different screening properties 
and dependences on frequency and temperature from the 
Meissner state, hence the results presented here are not 
valid for the mixed state. 

The paper is organized as follows. In SecHIl we give a 
brief introduction into vacuum fluctuation effects with 
emphasis on spin flip transitions and the Lamb shift. 
The expressions for the spin flip rate in isotropic and 
anisotropic planar multilayered structures is discussed 
in detail in Sec. Ill Al and in the Appendix. The results 
on atomic spin relaxation are given in Sec. IIIII with the 
screening properties investigated in Sec. IIII Al and the 
temperature dependence of the penetration depth investi- 
gated in Sec. HHBI Our conclusions are given in Sec. IIVI 



II. VACUUM FLUCTUATIONS 

The spontaneous decay of an excited atom is one of 
the most widely studied effects of the ground-state fluc- 
tuations of the electromagnetic vacuum {23?]. However, 



the presence of a dielectric or metallic body strongly 
modifies the vacuum field and its statistical properties 
24| , leading to a modification of the spontaneous emis- 
sion rate as well as the Lamb shift. In particular, the 
latter is associated with the appearance of an attrac- 
tive dispersion force, the Casimir-Polder force. The en- 
ergy shift is position-dependent and takes the form of 
Ai; = A(")£; + U{rA), where the first term A(")£: is 
the contribution to the Lamb shift in free space, while 
the second term U{rA) is the position-dependent van der 
Waals potential of an atom at position r^. Another re- 
lated issue associated to the modification of the vacuum 
field is the body-induced modification of the atomic spin 
fiip lifetime [101, which is the main subject of this paper. 



Due to the Meissner effect, the presence of a supercon- 
ducting material provide a shield for the field fluctuations 
arising from a conducting substrate. The energy gap of 



a superconductor is A„ 



ksTc, with ks the 



Boltzmann constant, and corresponds to the minimum 
excitation energy required to break a Cooper pair. In s- 
wave superconductors, the gap frequency is almost of the 
order of THz (e.g. ojg « 700 GHz for niobium [251), while 
in d-wave superconductors it is usually an order of mag- 
nitude larger due to the higher transition temperature 
(e.g. the BSCCO gap frequency is given in the literature 
as 7.5 THz [1^). All the physical phenomena resulting 
from field fluctuations are going to be attenuated as long 
as the range of frequencies involved are smaller than the 
gap frequency. 



The atomic spin flip lifetime is enhanced by the pres- 
ence of a superconductor as the typical transition fre- 
quencies are of the order of MHz and fall well below the 
gap frequency for both s-wave and d-wave superconduc- 
tors. The same cannot be said straightaway for the van 
der Waals potential. The induced energy shift and broad- 
ening results in the atom having a position-dependent 
polarizability associated with the relevant dipole transi- 
tions of the atom. The potential C/(r^) is given by the 
sum of two terms: a contribution resonant with the mag- 
netic dipole transitions U^{rA) and an off-resonant con- 
tribution t/°'"(r^) [23]. For an atom in its ground state, 
the resonant term vanishes and the off-resonant term is 
obtained by integrating over all the frequencies, which 
means that an attenuation of the potential for frequen- 
cies Lu < Ljg is not relevant on the overall effect. Hence, 
the Casimir-Polder force experienced by a ground state 
atom is largely unaffected by the presence of either s- 
wave or d-wave superconductors. However, the resonant 
term [/''(r^) may dominate f7°''(r^) for an atom in an 
excited state and a superconductor may play a signiflcant 
role. Typical optical transition frequencies for an atom 
are of the order of THz, and d-wave superconductors in 
principle may be able to alter their Casimir-Polder po- 
tential. 
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A. Spin flip rate 

The spin flip rate of an atom placed near a conducting 
or superconducting body is obtained within the formal- 
ism of macroscopic quantum electrodynamics as [lo| 



r = fiQ 



(3) 
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where (f\S\i) is the spin matrix element for the relevant 
transition, fis denotes the Bohr magneton, and 55 ~ 2 is 
the electron's g factor. The Green function G{r,r',uj) is 
the fundamental solution to the classical Helmholtz equa- 
tion and thus contains all necessary geometric and elec- 
tromagnetic information about the macroscopic bodies. 
We review in the Appendix how to evaluate the Green 
tensor for an anisotropic multilayered planar structure 
such as the one represented in Fig. [T] consisting of vac- 
uum, superconducting layer and metal substrate. 
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FIG. 1: Schematic representation of the geometry of the two 
structures considered: a three-layered geometry consisting of 
(a) vacuum-niobium-metal or (b) vacumm-BSCCO-metal. 

The spin flip rate for a planar multilayered structure 
made of isotropic layers is known to be [28| 
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where the spin matrix elements have been evaluated for 
the *^Rb ground state transition between the magnetic 
hyperfine sublevels |2,2) — > |2, 1). Let us recall here that 
Eq. Q is valid whenever the atomic transition wave- 
length can be regarded as the largest wavelength of the 
system. The generalized Fresnel reflection coefficient 
for a three-layer geometry reads 
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where d is the thickness of layer 2. The function r^f is 
the Fresnel reflection coefficient for TE waves at a planar 
interface defined as 
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with kf^ — kf — and kf — (a;^/c^)ei(aj) with the 
label i indicating the layer. The coefficient rjjf has the 
same form as ^fter the replacements ki^ k^z and 
k2z k^z have been made. 

The dielectric function e{uj) for a normal metal is given 
by e{uj) « 2ieo/fc^5^ where Sm is the skin depth of the 



metal and k = uj/c. In the two-fluid model To'], the di- 
electric function e{uj) of a superconductor can be written 
as 

1 . 2 



e(^) = 1 - 
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where S{T) = •\/2/a;/io(T„(r) and cr„(r) are the skin 
depth and the conductivity associated with the normally 
conducting electrons, respectively. The optical conduc- 
tivity corresponding to Eq. ([7]) is 
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Both the penetration depth and the skin depth are tem- 
perature dependent fioj, 
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where a is the electrical conductivity of the metal in the 
normally conducting state for T > Tc, and ns(T) and 
n„(r) are the electron densities of the superconducting 
and normal state, respectively, at a temperature T < 
Tc- The total electron density uq — ns{T) + nn{T) is 
constant, with na{T) — hq for T = and n„(T) — hq 
for T >Tc. It follows from Eq. (P that Eq. © can be 
written as 



A(0) 



\ nn{T)' 


-1/2 




-1/2 


no 









(10) 



1 for 



with a = 4 for an s-wave superconductor, and a 
a d-wave superconductor. 

In contrast, in a c?-wave superconductor the penetra- 
tion depth and thus the permittivity is highly anisotropic 
so that it can be written as a tensor e of the form 
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where tt and are the transverse and longitudinal scalar 
permittivities of the layer corresponding to A|| and A^, re- 
spectively. We calculate the spin flip rate for the BSCCO 
structure as 



ifJ-sgs)' 
8h 



dr^ l^[Sv'Bl}, + BU2kl] , (12) 



where B\}^ and Bj^2 are the relevant scattering coeffi- 
cients given in the Appendix. We assume that the whole 
structure of metallic substrate plus superconductor is in 
thermal equilibrium with its surroundings. The electro- 
magnetic field is then in a thermal state with temperature 
T, equal to the temperature of the materials. Therefore, 
both spin flip rates in Eq. Q and Eq. (fT^ need to be 
multiplied by a factor (rith + 1) where the mean thermal 
photon number is 



1 



„huj/kBT 
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III. ATOMIC SPIN RELAXATION 



superconducting layer as 



In this section, we study the atomic spin relaxation 
above a niobium and a BSCCO surface. We imagine that 
a *^Rb atom is held at a distance z from the supercon- 
ducting surface as shown in Fig. [TJ The spin flip lifetime 
T = 1/r is evaluated for the *^Rb ground-state transi- 
tion |2, 2) — |2, 1) with the transition frequency taken 
to be / = 560 kHz. The comparison between niobium 
and BSCCO is done in terms of the screening of the elec- 
tromagnetic field fluctuations affecting the atomic spin 
dynamics. Further comparison is done by studying the 
temperature dependence of the spin flip lifetime. 



A. Screening 

The alteration of the spin flip lifetime as a function 
of the relevant length scales of the system is an indica- 
tion of the capabilities of niobium and BSCCO to screen 
the electromagnetic field fluctuations that arise from the 
underlying metal substrate. We consider the multilay- 
ered structures as represented in Fig. [1] and we plot the 
atomic spin flip lifetime as a function of the distance z 
from the surface in Fig. [51 Two different thicknesses of 
the superconducting layer are taken into consideration. 
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FIG. 2: Atomic spin-flip lifetime near a niobium and a 
BSCCO surface plotted as a function of the distance from 
the surface 2 at T = 4.2 K. Continuos line: lifetime above 
niobium. Dashed line: lifetime above BSCCO. Two differ- 
ent thicknesses d of the superconducting layer are considered: 
d = 10 nm and d = \ jj.m. 

For both niobium and BSCCO, the lifetime increases 
with the thickness of the superconducting layer, and the 
longest lifetimes are obtained for the atom trapped above 
niobium (continuous line in Fig.[2|). During the supercon- 
ducting phase, magnetic field lines are expelled from the 
material except for a thin surface layer (with thickness 
of the order of the penetration depth), where shielding 
currents flow. Therefore, fluctuating magnetic fields are 
attenuated by the shielding currents. We now define a 
screening factor as a function of the thickness d of the 



where T{d = 0) is the spin flip lifetime associated with 
the bare copper substrate obtained by combining Eq. ^ 
and Eq. 

This screening factor S{d) is plotted in Fig. |3] for nio- 
bium and BSCCO films and for an atom held at a fixed 
10 fim distance from the surface. The screening S{d) in- 
creases in both structures for thicknesses d roughly the 
same size as the penetration depths at zero temperature 
(recall that A(0) = 35 nm for niobium and Ay (0) — 0.3 /im 
for BSCCO). For a superconducting layer that is thicker 
that the corresponding penetration depth, the screening 
effect saturates which confirms the fact that only the pen- 
etration depth provides an active screening of the electro- 
magnetic field fluctuations as already seen for dielectrics 

mm. 
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FIG. 3: Screening factor S{d) as a function of the thickness 
d for an atom held at a distance of 10 /im from the supercon- 
ducting surface. 

We emphasize that the results presented in this para- 
graph are only indicative for the two superconducting 
structures. The superconducting properties change both 
with the layer thickness being less than the coherence 
length (i.e. 30 nm for niobium at 4K). 



B. Temperature dependence 

Let us now investigate the temperature dependence of 
the spin flip lifetime. We have chosen a different thickness 
for each of the superconductors: 1 fim for niobium and 
2.5 /im for BSSCO, corresponding to the thickness pro- 
viding the maximum screening (see Fig. [3]) . As shown in 
Eq. Il]) and Eq. ^ , the penetration depth A is a function 
of the temperature, which is reflected in the temperature 
dependence of the optical conductivity in Eq. ([8|. How- 
ever, the spin flip lifetime is proportional to the mean 
thermal photon number nthi oc l/(nth + 1); which may 
become dominant for large enough temperatures. 
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The comparison is carried out in two stages. First, in 
Fig. |4] we show the hfetime t as a function of tempera- 
ture, while in Fig. [5] we show r as a function of T/Tc, 
in order to study how the different power laws affect 
the lifetimes near the transition temperature. The in- 




8.3 20 40 60 80 90 



Temperature T (K) 

FIG. 4: Spin-flip lifetime for the niobium and the BSCCO as 
a function of the temperature T. In the inset, the spin-flip 
lifetimes are plotted on a smaller temperature range 0.2 — 8.3 
K where T = 4.2 K refers to the liquid helium temperature. 
The atom-surface distance is fixed at z = 10 fim. 

set in Fig. m shows the spin flip lifetime for temperatures 
T < Tc^ where both niobium and BSCCO are super- 
conducting. The lifetime near superconducting niobium 
is several orders of magnitude larger than near super- 
conducting BSSCO. In particular, at the liquid helium 
temperature T = 4.2 K, the two lifetimes are approxi- 
mately Ts = 10^" s and = 5 10^ s as shown in the inset 
of Fig. m At the liquid nitrogen temperature T — 77 K, 
only BSCCO is still in the superconducting phase giv- 
ing a lifetime of roughly = 95 s, while the lifetime 
above (normally conducting) niobium is estimated to be 

Ts W 10-2 s. 
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FIG. 5: Spin flip lifetime for the niobium and BSCCO as 
a function of T/Tc- The range chosen in around T/Tc = 
1 in order to highlight the power law of the lifetime at the 
transition temperature. The atom-surface distance is fixed at 
z = 10 pm. 



In Fig. [SI the spin flip lifetime is plotted as a function 
of T/Tc based on Eqs. (|4]) and (fT2|) . For temperatures 
T/Tc > 1, we neglect the A dependence of the permit- 
tivity in Eq. O as n^(T) = 0. For T/Tc 1, the 
lifetime above BSCCO decreases with a smaller power 
law compared to niobium as a results of the penetration 
depth dependence on temperature. There is a marked 
difference in the transition from the normal state to the 
superconducting state, and we expect this to be visible 
in experiments comparing the two structures considered 
here. According to our results, the temperature depen- 
dence of the permittivity in d-wave and s-wave super- 
conductor can then be tested by observing the atomic 
relaxation rate near the transition temperature. 

We would like to draw the attention to the orientation 
of the atomic spin and how the spin flip rate changes 
according to it. The results presented in this section 
have been obtained for a randomly oriented spin. Due 
to the anisotropic permittivity of BSCCO, the shielding 
of the electromagnetic field fluctuations can in principle 
be different depending on which component couples to 
the atomic spin. However, our estimates of the spin flip 
lifetime for a spin oriented either perpendicular or par- 
allel to the interface suggest that the difference is not 
very pronounced. In fact, the spin flip lifetimes for the 
two orientations differ only for a factor of 2, the paral- 
lel orientation giving the longer lifetime, and it does not 
depend on the anisotropic permittivity. This indicates 
that, despite the transverse penetration depth being of 
the same order of magnitude as in metals, electromag- 
netic field fluctuations are shielded by screening currents 
independently from the anisotropy. 



IV. CONCLUSION 

We have considered the shielding of the electromag- 
netic field fluctuations provided by a superconductor. 
The modification of the vacuum field and its properties in 
the vicinity of a dielectric leads to the appearance of the 
Casimir-Polder potential and thermally-induced spin flip 
transitions. The frequencies associated with these two 
phenomena are very different and only when they fall 
below the superconducting gap frequency, the presence 
of the superconducting material is detectable. While the 
frequencies associated with spin transitions are always 
smaller than the gap frequency, the same can not be said 
for the frequencies involved in the Casimir-Polder poten- 
tial and numerical calculations need to be done. 

In this paper we have focused on the atomic spin re- 
laxation above a niobium and a BSCCO surface. We 
have chosen these two materials as representatives of an 
s-wave and a (high-temperature) d-wave superconduc- 
tor, respectively. We have considered a planar multilay- 
ered structure with the respective superconducting layer 
placed above a copper substrate. The shielding proper- 
ties of the two superconductors have been investigated 
by deflning a screening factor for the atomic spin relax- 
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ation and studying its dependence on the thickness of 
the superconducting layer. Despite not having obtained 
an analytical solution, we can still affirm that there is 
an active screening in the penetration depth layer of the 
superconductors. When the thickness is bigger than the 
penetration depth, the screening saturates. We believe 
that these results represent an important step towards 
the adoption of neutral atoms as sensitive probes in the 
study of properties of different types of superconductors. 

The temperature dependence of the permittivity in s- 
wave and d-wave superconductor follows two different 
power laws, such difference has been studied by com- 
paring the spin flip lifetime. When both niobium and 
BSCCO are in the superconducting phase, the niobium 
gives the longest lifetimes. The comparison has been car- 
ried on by investigating the spin flip lifetime near the su- 
perconducting transition. There is a marked difference 
between the two superconductors for T < Tc, which we 
believe can be observed experimentally to test the differ- 
ent temperature dependences. 
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with h = yjk"^ — if- and k = lo/c. The scattering dyadic 
Green function has a form similar to the unbounded 
Green function and can be formulated as 
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where in the lih. layer hij = ykf^ — with 
7 2 2 



(A4) 
(A5) 



and ft.11 — hi2 — h. The cylindrical vector wave functions 
are defined as 



Appendix A: Green function for anisotropic planar 
multilayers 
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In this section, we present our calculation of the Green 
tensor for an anisotropic multilayered structure based on 
the work of Li et al. in Ref. [2^. The geometry of the 
problem is shown in Fig. [TJd where the second layer is 
characterized by a tensor dielectric permittivity as in 
Eq. (jlip . The Green function will be expressed in terms 
of an expansion of the cylindrical vector wave functions 



adopting the cylindric basis {e^ , , 6;, 
cated at r' in the first layer (vacuum) 
can be written as 



}. For an atom lo- 
, the Green tensor 



G(r, r', uj) = Go(r, r', w) + G«(r, r', c.) , 



(Al) 



where Go(r, r',a;) is the unbounded (bulk) Green tensor 
representing the contribution of direct waves from the 
source at r' to the point r, while Gs(r,r',a;) describes 
the reflection of waves from the interface between the 
first and second layer. 



in which J{rir) is the Bessel function of the flrst kind 
and the first order, and fc,, is generally defined as fc,, = 

(?72 + /i2)l/2^ 

The coefficients B\\^, B]^2 ^-i^d B]^2 &re determined by 
satisfying the boundary conditions for the Green function 
at the interface between layers. For the sake of simplicity, 
we are now going to choose a reference system such that 
the coordinate r can be safely set equal to and the 
Bessel functions can be expanded for ri] — > such that 
Jo{rri) w 1 and Jn{rr]) w {rri/2)"- /n\. 

The spin flip rate as in Eq. ^ has been evaluated 
throughout this paper for a randomly oriented spin. We 
consider the relaxation rate to be the sum of the contribu- 
tion given by the atomic spin parallel or perpendicular 
to the planar substrate. We calculate the double curl 
of the Green function (^ x G(r, r') x ^')mn such that 
m,n — r, z, with the r component coupling to the paral- 
lel orientation of the spin and the z components coupling 
to the spin perpendicularly oriented. The double curl of 
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the scattering Green tensor reads 
^ X G,(r,r') X ^' = 



(A8) 
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The scattering coefficients for a general argument h / are 
given by 



-2ihid 



B 



M,N 



± 



1 + iif '^iif •^e'2''2c( 



i?f + R. 



H,V 2ih2d 



(A9) 



with 



H 



R 



RJ 



hf+i - hf 
hf+i + hf ' 

/l/[(i^i-a;2)/ij^i+ti)2fe/ + i] _ ^ 

/t/ + l [(OJI — tiJ2 ) 'i^ +^^2 fc^] 
fc/[(ail-lx>2)fe^^.l+<^2fc^^J 

hf+i{{uii-ui2)h?f+ui2k'j] 



1 



(AlO) 
(All) 



where d is the height of the second layer, and uj\ and uj^ 

arc the weighting coefficient that in our case equal to 1. 
To obtain the specific scattering coefficients, hf has to 
be substituted by hf \ or hf2- 
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